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ABSTRACT 



We consider a model in which the narrow hne regions (NLRs) of Seyfert 
galaxies are photoionized "in situ" by fast (300 - 1,000 km s~^), radiative 
shock waves driven into the interstellar medium of the galaxy by radio jets 
from the active nucleus. Such shocks are powerful sources of soft X-rays. 
We compute the expected ratio of the count rates in the ROSAT PSPC and 
Einstein IPC detectors to the [OIII]A5007 flux as a function of shock velocity, 
and compare these ratios with observations of type 2 Seyferts. // most of the 
observed soft X-ray emission from these galaxies originates in the NLR and the 
absorbing hydrogen column is similar to that inferred from the reddening of the 
NLR, a photoionizing shock model with shock velocity ~ 400 - 500 km s"^ is 
compatible with the observed ratios. High angular resolution observations with 
AXAF are needed to isolate the X-ray emission of the NLR and measure its 
absorbing column, thus providing a more conclusive test. We also calculate the 
expected coronal iron hne emission from the shocks. For most Seyfert 2s, the 
[Fe X]A6374/H/3 ratio is a factor of 2 - 14 lower than the predictions of 300 - 
500 km s~^ shock models, suggesting that less hot gas is present than required 
by these models. 



Subject headings: galaxies: active - galaxies: ISM - galaxies: jets - galaxies: 
nuclei - galaxies: Seyfert - shock waves 
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1. Introduction 

The gas in the narrow hne regions (NLRs) of Seyfert galaxies is known to be 
photoionized (e.g. Koski 1978; Ferland & Osterbrock 1986; Osterbrock 1989) and the 
prevaihng view is that the ionizing photons originate in a compact source, perhaps the 
accretion disk. However, strong associations between the narrow emission-hne and radio 
continuum properties of Seyfert nuclei have been found over the last 20 years. There are 
strong correlations between radio power and both line luminosity (de Bruyn & Wilson 1978) 
and hne width (Wilson & Willis 1980). Early VLA maps showed that the radio sources 
represent collimated ejection from the compact nucleus and that they have similar spatial 
scales and orientations to the NLR. These results led Wilson & Willis (1980) to suggest 
that the nucleus ejects radio components which interact with ambient gas and replenish the 
high kinetic energy and ionization of the NLR. Simple models of the momentum transfer 
between jets and ambient gas (Blandford & Konigl 1979) show that ambient gas with the 
obscrvationally-inferred NLR mass can be accelerated to the observed velocities (~ 10^ km 
s~^) in the gas crossing time of the NLR (~ 10^~^ yrs) for reasonable efficiency factors 
relating jet and radio powers (Wilson 1981, 1982). Subsequent imaging (e.g. Haniff et 
al. 1988; Bower et al. 1995; Capetti et al. 1995; Falcke, Wilson & Simpson 1998) and 
spectroscopic (Whittle et al. 1988; Axon et al. 1998) observations have confirmed that 
the structure of the NLR in many Seyferts is dominated by compression of interstellar gas 
by the radio ejecta and that these ejecta are an important source of "stirring" of the gas. 
Theoretical descriptions of this interaction have involved expanding radio lobes (Pedlar, 
Dyson & Unger 1985), bow shocks driven by the radio jets (Taylor, Dyson & Axon 1992; 
Ferruit et al. 1997) and the role of the jet cocoon (Steffen et al. 1997). 

The power required to accelerate the gas is typically ~ iq^i-^s ^^.g g-i Qygp ^ 10^~^ 
yrs to give the ~ iQ^^-^s kinetic energy present in the observed clouds of the NLR 
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(more kinetic energy could be present in gas too tenuous to be visible in line emission). 
The power radiated in observed line emission is somewhat higher at ~ 10^1-44 ^^.g g-i^ 
which is ^ 10^~'^ times higher than the radio luminosity (e.g. Fig. 1 of Wilson, Ward & 
Haniff 1988). Recently, Bicknell et al. (1998) have argued that the narrow line emission 
in Seyfert galaxies is, indeed, powered entirely by the radio-emitting jets. The jets are 
considered to drive several hundred to a thousand km s"^ radiative shocks into interstellar 
gas on the hundred pc scale. Such "photoionizing shocks" are powerful sources of ionizing 
radiation and create photoionized precursors (Daltabuit & Cox 1972), the optical spectra of 
which are similar to Seyfert galaxy NLRs (Dopita & Sutherland 1995, 1996 [hereafter DSI]; 
Morse, Raymond & Wilson 1996). In this picture, the conversion of jet kinetic energy to 
radio emission is much less efficient in Seyfert galaxies than in radio galaxies and radio-loud 
quasars; the Seyfert jets are thermally dominated while jets in radio-loud AGN may be 
dominated by relativistic particles and magnetic fields. 

It is very important to decide whether the energy source that powers the NLR is "in 
situ" mechanical motion (jets) or the more conventional compact photoionizing source. 
This issue is closely related to the fundamental question of whether the putative accretion 
disk loses energy primarily in mechanical (as jets or winds) or radiative form. 

Shock velocities of 300 - 1,000 km s~^, required to account for the emission-line 
spectra of Seyferts (Dopita & Sutherland 1995), generate gas of temperature 10^~^K and 
copious soft X-rays (cf. Laor 1998). This production of soft X-rays by the NLR is an 
inevitable prediction of the photoionizing shock model and can be used to distinguish it 
from photoionization by a hidden Seyfert 1 nucleus. Indeed, soft X-ray emission from the 
NLR has been isolated in NGC 1068 (Wilson et al. 1992), NGC 2110 (Weaver et al. 1995) 
and NGC 4151 (Morse ct al. 1995). For type 2 Seyfert galaxies, in general, the observed 
Einstein IPC soft X-ray flux (from the entire galaxy) and the [OIII]A5007 flux (from the 
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NLR) are similar (Fig. 1). This similarity is not inconsistent with photoionizing shocks 
being the power source, since ~ 4% (20%) of the kinetic power of a 300 (500) km s~^ 
radiative shock is radiated in the Einstein band (0.2 ~ 4 keV), while ~ 2% is radiated as 
[OIII]A5007 (DSI; Bicknell, Dopita & O'Dea 1997, hereafter BDO), and the expected X-ray 
spectrum is soft and photoelectric absorption can substantially reduce the observed flux. 
In this letter, we examine the relationship between soft X-ray and [Olll]A5007 emission 
expected in the photoionizing shock model, and compare it with Einstein and ROSAT 
X-ray observations. We also argue that the shocks should produce strong coronal iron hne 
emission, and compare the predicted and observed strengths. 

2. Comparison of [OIII]A5007 and Soft X-ray Fluxes 

2.1. Method 

The calculations of DSI provide the luminosity radiated in the [Olll]A5007 line per unit 
area of shock from both the post-shock gas and the photoionized precursor as a function of 
shock velocity and magnetic parameter. The pre-shock density is taken to be low enough 
that coUisional deexcitation is unimportant. We have used the model-predicted [OIII]A5007 
luminosities from DSI (with the fluxes reduced by a factor of 2, as we suspect they may 
be too high by this factor based on comparison with the models discussed below and an 
apparent factor of 2 error in equation 3.3 of DSI) and BDO. For the post-shock gas, the 
[OIII]A5007 luminosity was taken to be the average of the predictions for the four magnetic 
parameters considered by DSI. The exact value of the magnetic parameter is not critical 
because for shock velocities above ~ 200 km s~^, the [OIII] emission of the precursor 
dominates that of the post-shock gas. The [OIII]A5007 luminosity would be less than 
predicted if any of the following effects are important: a) the precursor is density bounded 
or the pre-shock gas presents a covering factor of < 1; b) the post-shock gas cools in an 
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unstable fashion, so that the fragmented gas does not intercept all of the ionizing flux (e.g. 
Innes, Giddings & Falle 1987a, b); or c) the pre-shock density is sufficiently high that the 
level of OIII suffers coUisional deexcitation. Of these, effect b) is unlikely to affect our 
results significantly in view of the small contribution of the post-shock gas to [OIII]A5007 
for the shock velocities of interest. Since the coUisional deexcitation density of the 
level is 7.0 x 10^ cm~^ (Osterbrock 1989), it is very unlikely that effect c) applies to the 
precursor gas, though post-shock gas could be coUisionally deexcited in the inner, denser 
parts of the NLR. The possibility that effect a) is significant means that the predicted 
[OIII]A5007 luminosities are most conservatively treated as upper hmits. 

We have also calculated the luminosity of the radiation emitted from the post-shock 
gas in the bands of the ROSAT PSPC (0.1 - 2.4 keV) and Einstein IPC (0.2 - 4 keV) 
detectors using an updated version of the radiative shock wave code described in Raymond 
(1979). The atomic rates are basically the same as those of the current version of the 
Raymond & Smith (1977) X-ray code used for the X-ray spectrum predictions in DSL The 
shock code includes time-dependent ionization and photoionization, but these have little 
effect on broad-band X-ray count rates. To obtain the total emission (over 47r sterad), we 
multiphed the X-ray flux radiated upstream by a factor of 2. This factor may be too large 
for soft photons because of photoelectric absorption by the downstream gas. The predicted 
ratio, R, of X-ray count rate to [OIII]A5007 flux can then be obtained as a function of shock 
velocity and effective hydrogen column density. The code also predicts the intensities of 
the [Fe X] and [Fe XIV] hues discussed below, including excitation by electron and proton 
impact, and by cascades following LS permitted excitations. 
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2.2. Results 

Figs 1 and 2 show the Einstein IPC and ROSAT PSPC count rates against the 
[OIII]A5007 flux. The dashed hues in each panel are loci of constant R obtained from the 
model as described above after attenuation by an equivalent hydrogen column of N^;^ = 
2Q20.5 (,j^-2^ which is the average Galactic column density towards the Seyfert 2 galaxies 
plotted in the flgures (see below). The [OIII] fluxes have been reduced by 0.2 mag, which 
is the expected obscuration for a normal gas to dust ratio. The solid lines are the same 
but with the model X-ray fluxes attenuated by N^^ — 10^^^^ cm~^, which is the column 
density corresponding to the average obscuration {Ay = 1.8 mag, derived from the Balmer 
decrement assuming intrinsic case B recombination values) of the NLRs, assuming a normal 
gas to dust ratio. The (narrow) hydrogen lines in Seyfert 2s come from a similar region to 
[OIII]A5007, so the model-predicted [OIII]A5007 fluxes have been reduced by 1.12 x Ay = 
2.0 mag. Because the model [OIII]A5007 fluxes may be overestimates, as discussed above, 
the predicted values of R are best treated as lower limits, so the predicted locations in Figs 
1 and 2 for a given shock velocity lie to the upper left of the plotted line. 

Also plotted in Figs 1 and 2 are observed X-ray count rates and [Olll]A5007 fluxes 
for type 2 Seyfert galaxies. Type 1 Seyferts and the so-called "narrow line X-ray galaxies" 
have been omitted (with the exception of NGC 2110, see flgure captions) because the 
equivalent hydrogen columns to their compact nuclei are sufficiently low that the soft X-ray 
emission is dominated by the Seyfert 1 nucleus. The equivalent hydrogen column towards 
the "hidden" Seyfert 1 nuclei in Seyfert 2 galaxies averages ~ 3 x 10^^ cm~^ (Turner et 
al. 1997), so the compact source is heavily attenuated in the Einstein and ROSAT bands. 
At the spatial resolution of the IPC and PSPC detectors, the measured X-ray count rate 
generally includes the entire host galaxy, the NLR and any residual transmitted or scattered 
emission from the compact nucleus. Thus these observed X-ray fluxes represent flrm upper 
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limits to the X-ray flux from the NLR. The [OIII]A5007 emission, on the other hand, was 
measured through a smaU aperture and is a good measure of the NLR hne emission. Thus 
the plotted data represent upper limits to the actual value of R in the NLR. 

Also plotted in Fig. 2 is a point representing the total PSPC count rate (after 
absorption) and [OIII]A5007 flux density (after obscuration) for a 3" (0.7 pc) length along 
the shock front of the LMC supernova remnant N132D (Morse et al. 1996). The plotted 
PSPC X-ray count rate was obtained by integrating the curve of 0.1 - 2 keV flux density 
(obtained from the ROSAT HRI image and corrected for absorption) versus distance 
perpendicular to the shock given in Fig. 11 of Morse et al. (1996). The predicted absorbed 
PSPC flux density was then found using the spectral model of Hwang et al. (1993) and 
the equivalent hydrogen column of N^ = 6.2 x 10^° cm~^, which is within a factor of 2 of 
the column used for the models represented by the dashed lines. The curve of [OIII]A5007 
flux density versus distance perpendicular to the shock, also given in Fig. 11 of Morse et 
al. (1996), was similarly integrated and was reduced by an obscuration Ay — 0.34 mag 
(equivalent to Nh = 6.2 x 10^° cm~^ for a normal gas to dust ratio). Fig. 2 shows that 
the resulting value of R for N132D is about 50 times higher than the average of the Seyfert 
galaxy points. Morse et al. (1996) argue convincingly that the X-ray and [OIII]A5007 
emissions from this region of N132D originate from a photoionizing shock of velocity 790 
km s~^, with some contribution to the ionizing flux from lower velocity shocks. The 
absorption-corrected ratio F(0.1 - 2 keV)/F([OIII]A5007) for this 0.7 pc length of shock is 
~ 45, with an uncertainty of about 25%. 

Examination of Figs 1 and 2 shows that the model generally overpredicts the observed 
R values (i.e. the predicted X-ray flux is too high) for high shock velocities (500 - 1,000 

km s^^) when only the Galactic column (N// = 10^" ''' cm^^) is included. For this column 
density, a shock velocity of ~ 300 km s~^ better describes the observed R values. 
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For Nh — 10^^-^ cm~^, the model of a 500 km shock is comparable to most observed 
detections and upper hmits, while the X-ray emission of a 300 km s~^ shock is ~ 2 orders of 
magnitude (the exact number being very sensitive to Nh) below the observed data points. 
If most of the observed soft X-ray emission originates in the NLR and the hydrogen column 
is similar to that inferred towards the optical (narrow) line emission, a photoionizing shock 
model with shock velocity ~ 400 - 500 km is compatible with the observations. On the 
other hand, if the X-ray flux from the NLR is actually substantially lower than the observed 
total emission from the galaxy, a lower value of the shock velocity would be indicated or a 
classical photoionization model favored. Future high spatial resolution, medium spectral 
resolution observations with the ACIS on AXAF will spatially resolve the X-ray emission 
of the NLR in nearby Seyferts, and thus separate it from both the galaxy disk and the 
compact nucleus. These measurements will also allow the gas temperature and N^^ to be 
determined, providing a much more precise check of the model. 



3. Coronal Iron Line Emission 

As emphasised by DSI, the post-shock gas should display a rich coUisionally ionized 
UV spectrum. However, these lines are very sensitive to reddening. Another prediction 
of fast shock models is strong emission in the coronal forbidden hues of iron, such as [Fe 
X]A6374 and [Fe X1V]A5303. These lines are expected whenever the shocked gas is hot 
enough to coUisionally ionize Fe to these species. Table 1 shows both absolute fluxes and 
fluxes relative to H/3 in these two hues. The H/3 emission (HPs) is that from only the 
post-shock gas - emission from the precursor is not included in the listed ratio. The range 
of shock velocities consistent with both the optical line ratios observed in Seyferts (Dopita 
& Sutherland 1995) and the soft X-ray observations (Section 2) is 300 - 500 km s~^. For 
such shock velocities, the ratio of total H/3 flux (UPr - shock plus precursor) to shock-only 
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flux (HPs) is 2.1 - 2.2 (DSI, tables 8 and 10). Thus the predicted [Fe X]A6374/H/?r ratios 
are in the range 0.1 - 0.19, excluding any [Fe X]A6374 from the photoionized precursor. 
The observed, reddening-corrected [Fe X]A6374/H/3 ratios in Seyfert 2s are in the range < 
0.014 - 0.16 (Koski 1978; Penston et al. 1984), with most galaxies in the range < 0.014 
- 0.06. The strongest [Fe X]A6374 emitters are consistent with a 300 - 500 km s~^ shock 
model, but most galaxies show an [Fe X]A6374/H/9 ratio which is a factor of 2 - 14 lower 
than the predictions. A similar conclusion comes from the general lack of detection of 
[Fe XIV]A5303 in Seyfert galaxies. Of course, a contribution from lower velocity (< 300 
km s"-*^) shocks would enhance the H/3 flux without significantly increasing [Fe X]A6374 
or [Fe XIV]A5303 and thus reduce the predicted ratio. However, the optical emission line 
spectrum would then be of low excitation and would not match a Seyfert NLR spectrum. 
Thus the weakness of the iron coronal lines may favor a conventional photoionization model 
over the shock model. Depletion of iron onto grains would reduce the coronal line fluxes, 
but it would decrease the X-ray emissivity as well, and much of the iron is liberated by the 
time radiative cooling sets in (Vancura et al. 1994). Thus the [Fe X]/H/3 ratios would be 
reduced by a modest amount in models with grains. 

This research was supported by NASA through grants NAG 53393 and NAG 81027, 
and by NSF through grant AST9527289. We are grateful to P. Ferruit, J. A. Morse and N. 
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Fig. 1. — A plot of the logarithm of the total Einstein IPC count rate (in counts s~^, left 
axis) versus the logarithm of the flux in [OIII]A5007 (in erg cm~^ ^~^)- The right axis gives 
the corresponding logarithm of the soft X-ray flux (in erg cm~^ s~^ in the 0.2 - 4.0 keV band, 
with no correction for absorption) . The points and upper limits are data for Seyfert 2 galaxies 
from Mulchaey et al. (1994). With the exception of NGC 2110, so called "narrow hne X-ray 
galaxies" are omitted. The vertical line descending from the upper NGC 2110 point joins 
the total observed count rate to the approximate count rate of extended X-ray emission 
some 4" N of the nucleus and within the NLR (from Weaver et al. 1995). The vertical hne 
descending from the upper NGC 1068 point joins the observed count rate from the entire 
galaxy to the approximate count rate of the spatially resolved emission in the ROSAT HRI 
image (from Wilson et al. 1992). The diagonal straight hues are predictions of photoionizing 
shock models, using calculations of [OIII]A5007 flux from either DSI (modified as mentioned 
in the text and used for shock velocities < 500 km s^^) or BDO (for a shock velocity of 1,000 
km s~^), and calculations of the Einstein IPC count rate described in the text. The shock 
velocities (in km s"^) are indicated. The dashed hues show the model predicted Einstein 
count rates reduced by the effects of photoelectric absorption by an equivalent hydrogen 
column of Nh = 10^°'^ cm~^ (which is the mean Galactic column to the galaxies plotted), 
and the model predicted [OIII]A5007 fluxes reduced by an obscuration A5007 = 0.2 mag. 
The sohd hues are the predictions of the model for an absorbing column of Nh — 10^^-^ 
cm~^ (which is the mean column inferred from the Balmer decrement of the NLR) and an 
obscuration A5007 — 2.0 mag. 

Fig. 2. — As Fig. 1, but for the ROSAT PSPC count rate. The data plotted are the Seyfert 
2 galaxies in Rush et al. (1996), with [OIII]A5007 fluxes from Whittle (1992). The rectangle 
represents a 0.7 pc length along the shock front of the LMC SNR N132D (see text). 
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Table 1. Coronal Line Emission from the Shock 



Shock velocity (km s ^) 


[Fe X]A6374 


[Fe XIV]A5303 


[Fe X]A6374/H/35 


[Fe XIV]A5303/H/35 


200 


2.16x10-^ 


0.0 


0.0025 


0.0 


300 


8.10x10-6 


1.90x10"^ 


0.41 


0.01 


350 


1.19x10"^ 


9.38x10-6 


0.38 


0.30 


400 


1.41x10-^ 


3.85x10-5 


0.35 


0.98 


500 


1.64x10-^ 


7.87x10-5 


0.23 


1.09 


1,000 


3.34x10-5 


1.52x10-^ 


0.13 


0.62 



The numbers in columns 2 and 3 represent the absolute fluxes (in erg cm ^ s ^) in the coronal lines 
emitted per unit area of shock for a pre-shock density of 1 atom cm"^. 
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